NONLINEAR DYNAMICS OF THE 3D PENDULUM* 
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Abstract. A 3D pendulum consists of a rigid body, supported at a fixed pivot, with three rotational degrees of freedom. 
The pendulum is acted on by a gravitational force. Symmetry assumptions are shown to lead to the planar ID pendulum and 
to the spherical 2D pendulum models as special cases. The case where the rigid body is asymmetric and the center of mass 
is distinct from the pivot location leads to the 3D pendulum. Full and reduced 3D pendulum models are introduced and used 
to study important features of the nonlinear dynamics: conserved quantities, equilibria, invariant manifolds, local dynamics 
near equilibria and invariant manifolds, and the presence of chaotic motions. These results demonstrate the rich and complex 
dynamics of the 3D pendulum. 
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1. Introduction. Pendulum models have been a rich source of examples in nonlinear dynamics and in 
recent decades, in nonlinear control. The most common rigid pendulum model consists of a mass particle 
that is attached to one end of a massless, rigid link; the other end of the link is fixed to a pivot point that 
provides a rotational joint for the link and mass particle. If the link and mass particle are constrained to 
move within a fixed plane, the system is referred to as a planar ID pendulum. If the link and mass particle 
are unconstrained, the system is referred to as a spherical 2D pendulum. Planar and spherical pendulum 
models have been studied in [Iffl [T] . 

Numerous extensions of simple pendulum models have been proposed. These include various categories 
of elastic pendulum models and multi-body pendulum models. Interesting examples of multi-body pendulum 
models are: a pendulum on a cart, an acrobot, a pendubot, a pendulum actuated by a reaction wheel, the 
Furuta pendulum, and pendula consisting of multiple coupled bodies. 

Pendulum models are useful for both pedagogical and research reasons. They represent physical mech- 
anisms that can be viewed as simplified academic versions of mechanical systems that arise, for example, 
in robotics and spacecraft. In addition to their important role in illustrating the fundamental techniques of 
nonlinear dynamics, pendulum models have motivated new research directions and applications in nonlinear 
dynamics. 

This paper considers a new 3D pendulum model, first introduced in [16j . and analyzes its nonlinear 
dynamical properties. This model consists of a rigid body, supported at a fixed pivot point that has three 
rotational degrees of freedom; it is acted on by a uniform gravity force. Control and disturbance forces and 
moments are ignored in this development. 

This paper arose out of our continuing research on a laboratory facility, referred to as the Triaxial 
Attitude Control Testbed (TACT). The TACT has been constructed to provide a testbed for a variety 
of physical experiments on attitude dynamics and attitude control. The most important feature of the 
TACT design is that it is supported by a three-dimensional air bearing that serves as an ideal frictionless 
pivot, allowing nearly unrestricted three degrees of rotation. The TACT has been described in several prior 
conference publications [H [9] . Issues of nonlinear dynamics for the TACT have been treated in [9l [8j . The 
present paper is partly motivated by the realization that the TACT is, in fact, a physical implementation of 
a 3D pendulum. 

2. Description of the 3D Pendulum. A rigid 3D pendulum is a rigid body supported by a fixed, 
frictionless pivot, acted on by gravitational forces. The supporting pivot allows the pendulum three rotational 
degrees of freedom. Uniform, constant gravity is assumed. The terminology 3D pendulum refers to the fact 
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that the pendulum is a rigid body with three spatial dimensions and the pendulum has three rotational 
degrees of freedom. 

Two reference frames are introduced. An inertial reference frame has its origin at the pivot; the first two 
axes lie in the horizontal plane and the third axis is vertical in the direction of gravity. A reference frame 
fixed to the pendulum body is also introduced. The origin of this body-fixed frame is also located at the 
pivot. In the body fixed frame, the moment of inertia of the pendulum is constant. This moment of inertia 
can be computed using the parallel axis theorem from the traditional moment of inertia with respect to a 
translated frame whose origin is located at the center of mass of the pendulum. Since the origin of the body 
fixed frame is located at the pivot, principal axes with respect to this frame can be defined for which the 
moment of inertia is a diagonal matrix. Note that the center of mass of the 3D pendulum may or may not 
lie on one of the principal axes defined in this way. 

Rotation matrices can be used to describe the attitude of the rigid 3D pendulum. A rotation matrix maps 
a representation of vectors expressed in the body-fixed frame to a representation expressed in the inertial 
frame. Rotation matrices provide global representations of the attitude of the pendulum, which is why they 
are utilized here. Other attitude representations, such as exponential coordinates, quaternions, or Euler 
angles, can also be used following standard descriptions, but each of the representations has a disadvantage 
of introducing an ambiguity or singularity. In this paper, the configuration of the rigid pendulum is a rotation 
matrix R in the special orthogonal group 50(3) defined as 

SO{3) = {i? e M.^''^ : RK" = /3X3, det(i?) ^ 1}. 

The associated angular velocity, expressed in the body-fixed frame, is denoted by w G M^. The constant 
inertia matrix, in the body-fixed frame, is denoted by the symbol J. The constant body-fixed vector from the 
pivot to the center of mass of the pendulum is denoted by p. The symbol g denotes the constant acceleration 
due to gravity. 

Three categories of 3D pendulum models are subsequently introduced and studied. The "full" dynamics 
of the 3D pendulum are based on Euler's equations that include the gravity moment and the rotational 
kinematics, expressed in terms of the angular velocity and a rotation matrix; this model describes the 
dynamics that evolves on TS0{3). Since the gravity moment depends solely on the direction of gravity in 
the pendulum fixed frame, it is possible to obtain a reduced model expressed in terms of the angular velocity 
and a unit vector that defines the direction of gravity in the pendulum fixed frame; this model describes the 
dynamics that evolve on TSO{3)/S^, and corresponds to the case of Lagrange-Poincare reduction Since 
there is a symmetry action given by a rotation about the gravity direction, Lagrange-Routh reduction [15] 
leads to a reduced model that is restricted to a constant momentum surface and is expressed in terms of the 
unit vector that defines the direction of gravity in the pendulum fixed frame and its derivative; this model 
describes the dynamics that evolve on TS^. Each of these 3D pendulum models provides special insight into 
the nonlinear dynamics. We develop each of these models in this paper, and we investigate the features of 
the nonlinear dynamics, namely invariants, equilibria, and stability for each model. 

3. 3D Pendulum Dynamics on TSO{3). The dynamics of the 3D pendulum are given by the Euler 
equation that includes the moment due to gravity: 

Jdi — Jlu X lu + mgp x R^e^. (3.1) 

The rotational kinematics equations are 

R = RQ. (3.2) 

Equations (|3.ip and (|3.2p define the full dynamics of a rigid pendulum on the tangent bundle TS0{3). 
In the above equations, 63 — (0,0, 1)^ is the direction of gravity in the inertial frame, so that R^e^ is the 
direction of gravity in the pendulum-fixed frame. The cross product notation a x b for vectors a,b gM.^ is 



ax [a2&3 - 0362, a^bi - 0163, ai&2 - a2&i] = ab, 



(3.3) 
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where, the skew-symmetric matrix a is defined as 
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A special case occurs if the center of mass of the rigid pendulum is located at the pivot. In this case 
p = 0, so that (|3.ip is given by Euler's equations with no gravity terms included. In the context of the rigid 
3D pendulum, this is referred to as the balanced case. Since there is a large literature on Euler's equations 
with no gravity moment and the associated rotational kinematics, this case is not considered further in this 
paper. Rather, the focus of this paper is on the unbalanced case, where p 7^ 0. 

3.1. Invariants of the 3D Pendulum. There are two conserved quantities for the rigid 3D pendulum. 
First, the total energy, which is the sum of the rotational kinetic energy and the gravitational potential energy, 
is conserved. In addition, there is a symmetry corresponding to rotations about the gravity direction through 
the pivot. This symmetry leads to conservation of the component of angular momentum about the gravity 
direction. These two results are summarized as follows. 

Proposition 1. The total energy 

and the component of the angular momentum vector about the vertical axis through the pivot 

h ^ Lu'^JR'^es. 



are each constant along motions of the rigid 3D pendulum. 

Proof. The proof follows by showing that the time derivative of the total energy and the time derivative 
of the angular momentum component about the vertical axis are each identically zero. We use p.ip and 
(13. 2p to compute the derivatives, yielding 

E — uj'^Juj — mgp^iTez — mguj'^{p x K^e^) + mgp'^{uj x K^es) — 0, 

and similarly, 

h^di'^JR'^es+ij'^JR'^es, 
= {Ju X LuyiR^ea) + {mgp x R-^e^YiR^e^) - uj'^J{lu x R^ea), 
= {JloY{uj X R-'es) - J(w x ii^eg) = 0. 

□ 

Conservation of the angular momentum component about the vertical axis will be revisited in Section 
[5] in the context of Noether's theorem, which states that the momentum map associated with the rotational 
symmetry about the gravity direction is conserved. 

3.2. Equilibria of the 3D Pendulum. To further imderstand the dynamics of the 3D pendulum, we 
study its equilibria. Equating the RHS of (|3.ip and (|3.2p to zero yields 

Jujc X LOc + mgp X R^es — 0, (3.5) 
i?cSc = 0. (3.6) 

Since i?c £ 50(3) is non-singular, and^ : M.^ K'^^'^ is a linear injection, i?c'2c = if and only if = 0. 
Substituting Wo = in p.5p . we obtain 



p X i?^e3 = 0. 



(3.7) 
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Hence, 

-^^063 = (3.8) 

or 

R>,^-^y (3.9) 

An attitude i?e is an equilibrium attitude if and only if the direction of gravity resolved in the body-fixed 
frame, i?o ^3, is coUinear with the vector p. If R]:e3 is in the same direction as the vector p, then (i?c, 0) is a 
hanging equilibrium of the 3D pendulum; if R^e^ is in the opposite direction as the vector p, then (Rc,0) is 
an inverted equilibrium of the 3D pendulum. 

Thus, if Rc defines an equilibrium attitude for the 3D pendulum, then a rotation of the 3D pendulum 
about the gravity vector by an arbitrary angle is also an equilibrium. Consequently, in TSO{3) there are 
two disjoint equilibrium manifolds of the 3D pendulum. The manifold corresponding to the first case in 
the above description is referred to as the hanging equilibrium manifold, since the center of mass is below 
the pivot for each attitude in the manifold. The manifold corresponding to the second case in the above 
description is referred to as the inverted equilibrium manifold, since the center of mass is above the pivot 
for each attitude in the manifold. 

Following (|3.8p and p.9p and the discussion above, we define 

[RU^{ReSO{3):R^e3 = /-:}, (3.10) 

[R],^\ReSO{3):R^e3^--f-r}, (3.11) 

as the hanging attitude manifold and the inverted attitude manifold, respectively. 
From dSiH]) and (j^ . 

H 4 |(i?,o) e rs'0(3) : i? e [i?],,}, (3.12) 

is the manifold of hanging equilibria and 

I = |(i?, 0) e TSOiS) : R e (3.13) 

is the manifold of inverted equilibria, and these two equilibrium manifolds are clearly distinct. 

3.3. Local Analysis of the 3D Pendulum near an Equilibrium. Consider a perturbation of the 
initial conditions from a hanging equilibrium (Rc, 0) of the 3D pendulum, using a perturbation parameter e. 
Let R'^ (t) and uj^ (t) represent the perturbed solution, corresponding to initial conditions R'^ (0) — Rc exp e6& 
and uj^{0) = eduj, where SQ,Suj G are constant vectors. Note that if e = then, (i?°(0), cj°(0)) = (i?c,0) 
and hence 

(i?"(t),c.O(t))^(i?e,0) (3.14) 

for all time t G M, which simply corresponds to the unperturbed equilibrium solution. 

Consider the solution to the perturbed equations of motion for the 3D pendulum. This solution satisfies 

Jlo' — Juj^ -K uo^ + mgp x {R'^yes, (3.15) 
7?^ = R%t)^. (3.16) 

Next, we differentiate both sides with respect to e and substitute e = 0, yielding 

Jcj° = Jw° X cj" + Jlu" X uj"^ + mgp x (i?°)^e3, (3.17) 
= (3.18) 
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where the subscripts denote derivatives. Substituting i?" = i?c and lj" = from (|3.14p into p.l7p and p.lSp 
yields 

Jtj", = mgp X (i?°)^e3, (3.19) 

ill = ReUJ°. (3.20) 

Now we define perturbation variables ALo{t) = uj'^{t) and AQ{t) ^ R]:R'^{t). It can be shown that ALu{t) — 
Ae(t). Thus, (|3?T9l) and ([3?20| can be written as 

jAe - ^^^^^Ae = 0. (3.21) 

\\p\\ 

^ 2 

Note that (13.2ip represents a mechanical system with mass matrix J, stiffness matrix f — but no 

\\P\\_ 

damping. Since p ^ is a negative-semidefinite matrix with two negative eigenvalues and one zero eigenvalue, 
the stiffness matrix is positive-semidefinite with two positive eigenvalues and one zero eigenvalue. The zero 
eigenvalue corresponds to rotations about the vertical axis, for which gravity has no influence. To see this 
more explicitly, we next perform a transformation of variables. 

As p ^ is a rank 2, symmetric, negative-semidefinite matrix, it follows from [51 [3] that one can simul- 
taneously diagonalize J and p^. Thus, there exists a non-singular matrix M such that J ~ MM"^ and 

— -rr-^p ^ = MAM^, where A is a diagonal matrix. Let A ~ diagimgli^ "^5^2, 0), where li and I2 are positive. 

\\P\\ 

Define x = Af^AG. Then expressing x = (xi, 2:2, 3:3) G M^, equation p.2ip can be written as 

xi + mglixi — 0, (3.22) 
X2 + mgl2X2 — 0, (3.23) 
ig = 0. (3.24) 

Thus, the variable x^ represents a perturbation in attitude of the 3D pendulum that corresponds to a rotation 
about the vertical axis. 

Due to the presence of imaginary and zero eigenvalues of the linearized equations, no conclusion can 
be made about the stability of the hanging equilibrium or the hanging equilibrium manifold of the 3D 
pendulum. Indeed, the local structure of trajectories in an open neighborhood the equilibrium is that of 
a center manifold; there are no stable or unstable manifolds. We next show that the hanging equilibrium 
manifold of the 3D pendulum is Lyapunov stable. 

Proposition 2. Consider the 3D pendulum model described by equations \8.1\) and \3.i3[) . Then, the 
hanging equilibrium manifold H given by h3.12\) is stable in the sense of Lyapunov. 

Proof. Consider the following positive-semidefinite function on TSO{2>) 

V{R,uj) = i uj^ Juj + mg{\\p\\ - p^K^e^,). (3.25) 

Note that V{R,{)) = for all {R,uj) G H and V{R,ljj) > elsewhere. Furthermore, the derivative along a 
solution of (|3.ip and p.2p is given by 

V{R, ijj) = Lo^ Jlu — mgp'^ R^ej, 

— uj^ {Juj X a; -I- mgp x R^e^) — mg p^ {—ujR^ e^) , 

= mg uj^ {p X R'^ es) + p^ {uj X R'^ es) =0. 



Thus, V is negative-semidefinite on TSO{3). Also, every sublevel set of the function V is compact. Therefore, 
the hanging equilibrium manifold H is Lyapunov stable. □ 
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Similarly, one can linearize the 3D pendulum dynamics about an equilibrium in the inverted equilibrium 
manifold. Expressing this linearization in terms of {xi, X2, X3) G M.^^ as in (|3.22p - (|3.24p . it can be shown that 
the linearization of the 3D pendulum about an inverted equilibrium can be written as 

xi — mghxi — 0, (3.26) 
X2 - mgl2X2 = 0, (3.27) 
X3 = 0. (3.28) 

The linearization of the 3D pendulum about an inverted equilibrium results in a system that has two positive 
eigenvalues, two negative eigenvalues and two zero eigenvalues. Thus, the inverted equilibrium has a two 
dimensional stable manifold, a two dimensional unstable manifold and a two dimensional center manifold. 
It is clear that due to the presence of the two positive eigenvalues, the inverted equilibrium is unstable. 

Proposition 3. Consider the 3D pendulum model described by equations 13. 1\) and i3.2l) . Then, each 
equilibrium in the inverted equilibrium manifold I given by !i3.13]} is unstable. 

4. Lagrange Poincare Reduced 3D Pendulum Dynamics on TS0{3)/ . The equations of 
motion (|3.ip and (j3.2p for the 3D pendulum are viewed as a model for the dynamics on the tangent bundle 
TSO{3) [S]; these are referred to as the full equations of motion since they characterize the full attitude of 
the rigid pendulum. As there is a rotational symmetry corresponding to the group of rotations about the 
vertical axis through the pivot and the associated angular momentum component is conserved, it is possible 
to obtain a lower dimensional reduced model for the rigid pendulum. This Lagrange-Poincare reduction is 
based on the fact that the dynamics and kinematics equations can be written in terms of the reduced attitude 
vector r = R^e^ G S^, which is the unit vector that expresses the gravity direction in the body-fixed frame 

M- 

Specifically, denote the group action of 9 e on 50(3) by $e : 50(3) 50(3), ^eiR) = exp(6'e3)i?. 
This induces an equivalence class by identifying elements of 50(3) that belong to the same orbit; explicitly, 
for i?i,i?2 e 50(3), we write Ri ~ R2 if there exists a. 9 G such that $e(i?i) = i?2- The orbit space 
50(3)/5^ is the set of equivalence classes, 

[R] A {^g{R) e 50(3) : e e 5I}. (4.1) 

For this equivalence relation, it is easy to see that Ri ^ R2 if and only if RJes = R^c-a and hence the 
equivalence class in (14. ip can alternately, be expressed as 

[R] ^ {R, e 50(3) : i^Jeg = R^e^}. (4.2) 

Thus, for each R e 50(3), [R] can be identified with T = R'^es e and hence SO{3)/S^ = 5^. This group 
action induces a projection H : 50(3) -> 50(3)/5^ = 5^ given by U{R) = R^e^. 

Proposition 4 ( P?] ) . The dynamics of the 3D pendulum given by Ii3.1\) and i3.S^) induce a flow on the 
quotient space TSO{3)/ , through the projection tt : r50(3) TSO{3)/ defined as tt{R, fi) = {R^e^, fi), 
given by the dynamics 

Jlo — Juj X uj + mgp x T, (4.3) 
and the kinematics for the reduced attitude 

t = r XLj. (4.4) 

Furthermore, TS0{3)/S^ ^ 5^ x M^. 

Equations (|4.3p and (|4.4p are expressed in a non-canonical form; they are referred to as the reduced 
attitude dynamics of the 3D pendulum on T50(3)/5^. 



Nonlinear Dynamics of the 3D Pendulum 



7 



4.1. Special Cases of the 3D Pendulum. Three interesting special cases are now examined. Suppose 
that the 3D pendulum is axisymmctric, that is two of the principal moments of inertia of the pendulum are 
identical and the pivot is located on the axis of symmetry of the pendulum. Assume that the body-fixed axes 
are selected so that J — diag(Jt, Jt, Ja) and p — PsC^ where Ps is a positive scalar constant. Consequently, 
equation (|4.3p can be written in scalar form, as 

' JtLO^ = {Jt - Ja) UJyLOz - mgpsTy, 

< JtLOy = (Ja - Jt) ijJz'jJx + mgpsVx, (4.5) 

^Ja^^z^ 0. 

From the last equation in (|4.5p . we see that the component of the 3D pendulum angular velocity vector 
about its axis of symmetry is constant. This means that a constant ujz defines an invariant manifold for 
the 3D pendulum dynamics. The special case that uj^ = c, where c G R, leads to invariant dynamics of the 
axisymmetric 3D pendulum, described as follows: 

Proposition 5. Assume the 3D pendulum has a single axis of symmetry and the pivot and the center 
of mass are located on the axis of symmetry of the pendulum as above. The equations of motion of the 3D 
pendulum define an induced flow on the fiber bundle x corresponding to uJz — c, where c ^ 0, given by 
the equations 

Jt LOx = c{Jt - Ja) LOy " TUgPsTy, 

Jt uJy = c{Ja - Jt)ujx + mgpsT^, 

r = r X [ujx LOy ■ 

These equations describe the dynamics of a spinning top where c S M denotes the spin rate. 

If the spin rate of the 3D pendulum about its axis of symmetry is zero, the following result is obtained. 

Proposition 6. Assume the 3D pendulum has a single axis of symmetry and the pivot and the center 
of mass are located on the axis of symmetry of the pendulum as above. The equations of motion of the 3D 
pendulum define an induced flow on the fiber bundle x corresponding to tOz = 0, given by the equations 

JtUJx = -mgpsTy, 
JtLOy = mgpsT^, 

t^T x[uj, LOy 0]^. 

These equations describe the dynamics of a 2D spherical pendulum. 

Now assume the spin rate of the 3D pendulum about its axis of symmetry is zero and, in addition, 
LOx{^) = and ry(0) = for the axisymmetric 3D pendulum. Then Proposition [S] yields invariant dynamics 
of the axisymmetric 3D pendulum such that uJx{t) = ^y{t) = for all i > 0. Therefore T can be parameterized 
by an angle as F = [— sin 6* cos6']'^. This yields the following result. 

Proposition 7. Assume the 3D pendulum has a single axis of symmetry and the pivot is located on 
the axis of symmetry of the pendulum as above. The equations of motion of the rigid pendulum define an 
induced flow on the tangent bundle TS^ corresponding to lo^ = 0, lo^ ~ given by the equations 

Jtujy = -mgpsSinO, 

9 = LOy. 

These equations describe the dynamics of a ID planar pendulum. 

Thus for an axially symmetric 3D pendulum with the pivot located on the axis of symmetry, the well 
known 2D spherical pendulum and the ID planar pendulum can be viewed as special cases of the 3D 
pendulum dynamics. For an axially symmetric 3D pendulum with the pivot located on the axis of symmetry, 
the induced dynamics corresponding to a nonzero constant value of lOz is fundamentally different from the 
dynamics of the spherical pendulum; these dynamics seem not to have been previously studied. It should be 
emphasized that if the 3D pendulum is asymmetric then the dynamics are general in the sense that neither 
the 2D spherical pendulum dynamics nor the ID planar pendulum dynamics are special cases. 
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4.2. Invariants of the Lagrange— Poincare Reduced Model. In a previous section, we obtained 
two integrals of motion for the full model of the 3D pendulum. In this section we summarize similar results 
for the Lagrange-Poincare reduced model of the 3D pendulum considered here. 

Proposition 8. The total energy 

E=^uj'^Juj-mgp'^r, (4.6) 

and the component of the angular momentum vector about the vertical axis through the pivot 

h = uj'^JT. 

are each constant along trajectories of the Lagrange-Poincare reduced model of the 3D pendulum given by 
and CT . 

4.3. Equilibria of the Lagrange Poincare Reduced Model. We study the equilibria of the 
Lagrange-Poincare reduced equations of motion of the 3D pendulum given by (|4.3p and (|4.4p . Equating 
the RHS of (14. 3|) and ()4.4|) to zero yields that a natural equilibrium (FcWe) satisfies 

Jujc X cjc + mgp X To = 0, (4.7) 
r.XLj,^ 0. (4.8) 

Equation (|4.8p implies that uje — kT^, where fc e R, and substituting this into (|4.7p yields 

fc2 JEc X Pc + mgp X Pe 0. (4.9) 

Note that depending on whether k is equal to zero or not, one can obtain static {lOc = 0) or dynamic (lUc ^ 0) 
equilibrium. The dynamic case corresponds to relative equilibria of the 3D pendulum. 

Without loss of generality, we assume that the moment of inertia matrix is diagonal, i.e. J — diag( Ji, J2, J3), 
where Ji > J2 > J3 > 0. This is achieved by choosing the body- fixed reference frame such that the body- fixed 
axes lie along the principal axes. Note that the vector from the pivot point to the center of mass p may not 
lie along any principal axis. The following result describes the equilibria structure of the Lagrange-Poincare 
reduced equations. 

Proposition 9. Consider the Lagrange-Poincare model of the 3D pendulum given by j| ) and 
The equilibria (TcLUc) of the Lagrange-Poincare model are given as follows. 

1. The hanging equilibrium: ^jf^jO^ ^ 

2. The inverted equilibrium: ("jCT'O^j 



3. Two relative equilibria: 



^-v r^,-v).f-^,-,/^.-vu. (4.0) 



4- One dimensional relative equilibrium manifolds: 



sgn(a)^,./^n„ ) , (4.11) 



where — {J — -^Isxa)^^ P G 1^'^ corresponding to a G Ci, i G {1,2,3,4} and 

A-(-^,o)u(J^,^), £2^(0,^^), c,^ij^,j-^), c^^ij^,j-^). 

Also, sgn(-) denotes the sign function. 

One dimensional relative equilibrium manifolds: 



\\na\ 



-sgn(a)^,-.//^n„), (4.12) 



where Ua and a are as given above in Case 4- 
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The families of equilibria given in ^J^Tl^ and {4-12^ converge to the hanging equilibrium and the inverted 
equilibrium when a —>■ 0, and they converge to the third equilibria given in when a ±oo. // the 

vector from the pivot to the center of mass p, lies on a principal axis, i.e. p x Ci — for some i G {1, 2, 3}, 
then 1^. j j[ ) and can be rewritten as {(6^,76^), (— ej,7ej)} for 7 g R. 

Furthermore, there exist additional equilibria under the following assumptions on the moment of inertia 
matrix J and the vector from the pivot to the center of mass p — [pi p2 Ps]^ . 

6. Ji, J2 and J3 are distinct and pi — 0. Then there exist one- dimensional relative equilibrium mani- 
folds: 



(4.13) 



lb.||'Vlb«ll V'V b.ll' VIIp' 

for^^ {1,2,3}, where p, = (7,^,^), p, = (^,7,^), p, ^ (^,^,^) and 

7 e M. 
7. Ji = J2 7^ Js- 

(a) If pi — P2 — O. Then there exist two-dimensional relative equilibrium manifolds: 

g j mg \ / q j mg 

M'iMV'K Ikll' Vlkll'^, 

where q = (7, 6, jflj^ ) and 7, (5 G K. 

(b) If P3 — 0. Then there exist one- dimensional relative equilibrium manifolds: 

r I mg \ / r I mg 

where r = (7^^, -jfrj^.l) and 7 £ M. 
Proof. From (|4.9p . an equilibrium (FcCJo) satisfies 

/cVFo + m^p = fciFe, (4.14) 

for a constant fci e R. We solve this equation to obtain the expression for an equilibrium attitude Fe for 
two cases; when fci = and when fci ^ 0. The corresponding value of the constant k yields the expression 
for the equilibrium angular velocity as lOc = kVf,. 

Equilibria 3: Suppose fci = 0. It follows that fc ^ from (|4.14p . Thus, we have Fc = -^J^^p. Since 
llFoll — 1, we obtain = mg\\J^^ p\\, which gives (|4.10p . 

Equilibria 1, 2, 4 and 5: Suppose fci 7^ 0. If fc = 0, (|4.14p yields the hanging and the inverted equilibrium. 
Suppose k^O. Define a £ IR\{0}, and v = fciFc £ R'\ Then, can be written as 

(a J - I3x3)v = -rngp. (4-15) 

Note that for a G R\{0, j^, -j^, -j^} the matrix (J— ^/sxs) is invertible. Then, (|4.15p can be solved to obtain 
V — —^{J — ^13x3)^^ P- Since ||Fc|| = 1, we have ||w|| = ||fciFc|| = |fci|. We consider two sub-cases; when 
fci > 0, and fci < 0. 

If fci > 0, we have fci — \\v\\ and a > 0. Thus, we obtain the expression for equilibria attitudes as 

= TTir " "TT^ " ~^g^(")Tr^- (4.16) 

M IM \\na\\ 

where ^ [J — ^^3x3)^^P G ■ Since fc^ = afci = a \\v\\, we obtain the expression for equilibria angular 
velocities as 



= fcFe = TVaK' n ||=TW|| 4.17) 
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Thus, (|4.16p and (|4.17p correspond to the families of equihbria given by (|4.1ip and (|4.12p for a > 0. Consider 
the limiting case when a — > cx3. We have 



lim 



lim 



(J-/3x3/a)-V 



^3x3 



/a)-V|| 



Similarly, 



lim 

a — >oo 



mg 



mg 



\J-'P\\ 



Thus, as a — > oo, (|4.16p and (|4.17p converges to the first relative equilibria given in (|4.10p . 

Similarly, if fci < 0, we have fci = — a < 0, and k'^ — — q;||u||. Thus, the relative equilibria are 
described by 



V 

M 



WW 



-Sgn(Q;) 



WW 



± 



mg 



(4.18) 



which corresponds to the families of equilibria given by (|4.1ip and (|4.12p for a < 0. It can be similarly 
shown that they converges to the third relative equilibria given by (j4.10p as a ^ — oo. 

Next, consider (|4.1ip and (|4.12p as a — > 0. Expressing Ua = a{aJ — l3x3)^^P, we obtain 



lim 



-sgn(Q!) 



\\na\ 



lim 

Q-+0 + 



which corresponds to the inverted attitude. Similarly, 



lim 

q;-»0- 



-sgn(a)- 



lim 



|M-/3X3)-VII 



(aJ -/3X3) V 



a^O- ||(aJ-/3x3)"VII 

which corresponds to the hanging attitude. For the angular velocity term, 

P 



P 

\\P\\ 



_P_ 

\\P\\ 



lim 



mg 

¥\\' 



lim ^/mg\\rl\\ 



n 

R 



i^T—rr lim -v/ mgllapll = 0. 

p\\ a^O 



Thus, as a —>■ 0~, (14. lip and (|4.12p yields the hanging equilibrium given in case (1), and similarly, as a 
(|4.1ip and (|4.12p yields the inverted equilibrium given in case (2). 

Now suppose the vector from the pivot to the center of mass p, lies on a principal axis, i.e 
Bi = for some i G {1,2,3}. Then the vector p can be expressed as p = se^, where s £ M. 
for all a G M\{0, -j^ , , -j^}, (aJ — I^xs) is an invertible diagonal matrix, and hence Ua = ^{c^J 



0^ 



p X 
Then 



3x3) 



P 



Then, it follows from (|4.1ip and (|4.12p that the equilibria can be rewritten 



a{aJi — 1) 
as {(e,,7e^), (-e,,7e i)} fo r 7 £ R. 

The solution of (|4.15p for a e M yields all possible equilibria of (|4.3p and (|4.4p . Equation (|4.1ip and 
(|4.12p present the solution of (|4.15p for all a G M\{0, -j^, -^}. As shown before, a = yields the hanging 

and the inverted equilibrium. If a G {-j^, 'J^jjz}^ (|4.15p can have solutions under certain conditions. This 
yields the additional equilibria of (|4.3p and (|4!4p given in Case 6 and 7. 

Equilibria 6: Suppose Ji, J2 and J3 are distinct. Then it is easy to see that for a = 1/Ji, (|4.15p has a 
solution iff pi = 0. In this case, (|4.15p can be written as 














"0" 





J2-J1 





V — —mgJi 


P2 








J3-J1 




P3 



Since a > 0, it can be shown as in (|4.16p that To 



Pi 



and ujc — i 



llPill V Ibii 

and 7 G M. Similarly, one can yield solutions of (|4.15p for the case where a = I/J2 and a = I/J3 iff p2 = 



y^pi, wherepi = (7, j^, j^) 
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and /93 = 0, respectively. Thus for distinct principal moments of inertia, one obtains the equilibria given in 
case 5, corresponding to the condition pi — where i G {1,2,3}. 

Equilibria 7: (a) Suppose Ji = J2 ^ Js- Then it is easy to see that for a ~ 1/Ji = I/J2, (|4.15p has a 
solution iff pi = P2 = 0. In this case, (|4.15p can be written as 














"0" 











V = —mgJi 











h-Jl 




P3 



Since a > 0, it can be shown as in (I4.16P that Fc 



and 



-q, where q = (7, S, 7-377) and 



7,(5 e M. 

(b) Similar to above, one can yield solutions of (|4.15p for the case where a = I/J3 iff = 0. 

Thus for Ji = J2 ^ J^, one obtains the equilibria given in case 7(a) and 7(b) under the specified 
conditions. Finally for the case where Ji = J2 = J3, there are no additional solutions of (|4.15p for a G 
I— — — T □ 

Numerical example. We show the equilibrium structure of a particular 3D pendulum model. We choose 
an elliptic cylinder with its semimajor axis a = 0.8 m, semi-minor axis b = 0.2 to, and height 0.6 m. The 
pivot point is located at the surface of the upper ellipse, and it is offset from the center by [— |, 0]. The 
moment of inertia is given by J = diag(0.4486, 0.3943, 0.0772) and the vector from the pivot to the mass 
center is p = [-0.0140, 0.1044, 0.4989]. One of hanging attitudes is shown in Fig. |4.1(a) 

Figures 4.1(b) 4.1(d) show the relative equilibria attitudes on S^, where the top corresponds F = —63, 



and the bottom corresponds to F = 63. The inverted equilibrium is denoted by a red dot, and the hanging 
equilibrium is denoted by a blue dot, and the equilibria of (|4.10p are located at the intersection of the blue 
line and the red line in Fig. 4.1(c) The families of the relative equilibria given by (|4.1ip and (|4.12p are shown 



by the four segments of solid lines corresponding to £j, where i G {1, 2, 3, 4} and the value of a is represented 
by color-shading; a varies from —00 (blue color) to 00 (red color). Note that the reduced attitude in both 
(|4.1ip and (I4.12p . are the same and these families of equilibria only differ by a sign in the angular velocity 
vector at the equilibrium. 

The relative equilibria for a G (—00, 0) are shown by a segment of a blue line in Fig. 



4.1 c 



which starts 

from the third equilibria given by (|4.10p . and converges to the inverted equilibrium. For a G {0,-j^), three 



disjoint segments of the relative equilibria attitudes are shown in Fig. 4.1(d) a G (0, jj-) — C2 (upper- left 



moving counter-clockwise from blue to cyan), a G {j^j j^) ^ ^3 (upper-right, moving counter-clockwise from 



cyan to green), and a G (771 7^) = ^4 (lower center, moving upward from green to orange). The relative 



equilibria attitudes for a G {-j^,oo) are shown in Fig. 4.1(c) by a red line segment, which converges to the 
blue line at the third equilibrium given by (j4.10p . 

Since no component of the center of mass vector vanishes, there are no additional equilibria. In summary, 
the hanging attitude, the inverted attitude, and the attitude given by (j4.10p are equilibrium attitudes, and 
there are four mutually disjoint equilibrium attitude segments corresponding to Ci, where i G {1,2,3,4}. 

Relation to the equilibrium manifolds of the full model. Let 

Then it follows from Proposition [5] that (F/uO) and (F^jO) are equilibria of the Lagrange-Poincare reduced 
model of the 3D pendulum. These are called the hanging equilibrium and the inverted equilibrium of the 
Lagrange-Poincare reduced model, respectively. 

Let (i?e, 0) denote an equilibrium in either the hanging equilibrium manifold or the inverted equilibrium 
manifold of the full equations (|3.ip and (|3.2p and tt : TSO{i) — > TSO{3)/ be the projection as in 
Proposition m Then, it can be shown that either 7r(_Rc,0) — (F/i,0) or 7r(i?c,0) = (Fi,0). Thus, the hanging 
and the inverted equilibrium manifold of the full equations are identified with the hanging and the inverted 
equilibrium of the Lagrange-Poincare reduced equations. 
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(a) A hanging equilibrium (b) 3D view (c) Top view (d) Bottom view 

Fig. 4.1. Relative equilibria attitudes for an elliptic cylinder 3D pendulum model 

To study the properties of the equihbrium manifolds, it is advantageous to consider the hanging and the 
inverted equihbrium of the equations of the 3D penduhim in terms of the reduced attitude as in (j4.3p and 
(|4.4p . The following result provides the identification. 

Proposition 10 ([?]). The hanging equilibrium manifold and the inverted equilibrium manifold of the 
3D pendulum given by i3.1\) and i3.2\) are identified with the hanging equilibrium {Tfi,0) and the inverted 
equilibrium (r^jO) of the reduced attitude equations given by j|) and |^..^[). 



4.4. Local Analysis of the Lagrange Poincare Reduced Model near an Equilibrium. In the 

last section, we showed that the Lagrange-Poincare reduced model of the 3D pendulum can have equilibria 
with non-zero angular velocities. Also, there are only two static equilibria, namely the hanging equilib- 
rium and the inverted equilibrium. As stated in Proposition llOi these equilibria correspond to the disjoint 
equilibrium manifolds of the full equations of the 3D pendulum. 

We next focus on these static equilibria of the Lagrange-Poincare reduced equations. The identification 
mentioned in PropositiontTOlrelates properties of the equilibrium manifolds of the full equations and the equi- 
libria of the Lagrange-Poincare reduced equations. We deduce the stability of the hanging and the inverted 
equilibrium manifolds of the full equations by studying the stability property of the hanging equilibrium and 
the inverted equilibrium of the Lagrange-Poincare reduced equations. 

Consider the linearization of (I4.3I) - (I4.4I) about an equilibrium (r/i,0) = {H^e^jO), where (i?o,0) is an 

— 5, the linearization of (14.31)- 



equilibrium of the hanging equilibrium manifold H. Since dim TSO{3)/ 

dMl) about (Th,0) evolves on R^. 

Proposition 11. The linearization of the Lagrange-Poincare reduced equations for the 3D pendu- 
lum, about the equilibrium {Th,0) = {R^e3,0) described by equations can be expressed using 
(a;i,a;2,ii,a;2,a;3) e according to 13.2^) - (J7SJ\ I. 

Proof. Consider a perturbation in terms of the perturbation parameter £ S R as before. Let oj'^{t)) 
denote the perturbed solution of l|3.ip - (|3.2[) . Since F — R^es, the perturbed solution of (|4.3|) - (|4.4[) is 
given by {T^ (t) , oj'^ (t)) where r^(t) — {R'^y{t)e3. Define the perturbation variables Aoj{t) = uj'^it) and 
AT{t) = T°{t) = {R°y{t)e3. From definition of AO in Subsection [331 note that 

AF = -AGiJjea = F^Ae e Tr, S^. 

Then from p.l9|) and the definition of AF, it can be easily shown that the linearization of (I4.3p - ()4.4|) is 
given by 

JAuj ^mgp AT, (4.19) 
AF = f,, Aw. (4.20) 

Next, we express (j4.19p and (|4.20p in terms of {x,i). Specifically, we show that (AF, Alj) G Ir^S'^ x R-^ can 
be expressed using (xi, a;2, ii, ±2, 2:3) G R'^. 
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Since x = M'^AQ, and M is nonsingular, Auj — M~''^x and AF — ThM~'^x. We now give an orthogonal 
decomposition of the vector AO — M~'^x into a component along the vector p and a component normal to 
the vector p. This decomposition is 



M"-x = -^{M-^x) + ^ p-{M--x) 

\p\r WpW 



where 



p 2 ^ 

p'^{M "^x) p e span{p} and — (M "^x) £ span{p}^. Thus, 



ppr ^ - Il^ll 

Ar = Ae = = i_pMAx, 

\\P\\ rng\\p\\^ 

does not depend on ^3 since A = diag(m(//i, to(7Z2, 0). Thus, we can express the linearization of 
at (r;i,0) = (i?Je3,0) in terms of the variables {xi,X2,ii,X2,X3) according to (|3.22p - (|3.24p . □ 

Remark 1. Note that due to our careful choice of variables, one can discard X3 from (|3.22p - (|3.24[) when 
studying the stability properties of the inverted equilibrium manifold. Thus, X3 corresponds to a component 
of the perturbation in the attitude that is tangential to the inverted equilibrium manifold. However, the 
angular velocity corresponding to X3 given by is is retained. 

Summarizing the above, the linearization of l|4.3p - l|4.4p about the hanging equilibrium (F^ , 0) is expressed 

as 

Xi + mghxi = 0, (4.21) 
X2 + mgl2X2 = 0, (4.22) 
X3 = 0. (4.23) 

It is clear that due to the presence of zero and imaginary eigenvalues, one cannot arrive at a conclusion 
about the stability of the hanging equilibrium (r;i,0) from the linear analysis. Therefore we next consider 
Lyapunov analysis. 

Proposition 12. The hanging equilibrium {Th,0) — (n^'O^^ '^f reduced dynamics of the 3D 
pendulum described by equations li4-3\ ) and is stable in the sense of Lyapunov. 

Proof. Consider the candidate Lyapunov function 

V{r, Lu) = ^ uj^Jlu + mg{\\p\\ - p^r). (4.24) 

Note that V{rh, 0) = and V{T, a;) > elsewhere. Furthermore, the derivative along a solution of (|4.3p and 
(|4.4p is given by 

V^(F, Lu) — lJ^ Juj — mgp"^T, 

= lJ^ {Jlo X ciJ + mgp x F) — ragp"^ {T x w), 
= ui^mgp X F — mgp^T x cj = 0. 

Thus, the hanging equilibrium is Lyapunov stable. □ 
Remark 2. Note that combining Proposition[T2]with Proposition [10] immediately confirms the stability 
result for the hanging equilibrium manifold in Proposition [2j 

We next study the local properties of the Lagrange-Poincare reduced equations of the 3D pendu- 
lum near the inverted equilibrium (Fj,0). Consider the linearization of ()4.3p - (l4.4|) about an equilibrium 
(ri,0) = (7^363,0), where (i?o,0) is an equilibrium of the inverted equilibrium manifold I. A result similar 
to Proposition [TT] follows. 

Proposition 13. The linearization of the Lagrange-Poincare reduced equations for the 3D pendu- 
lum, about the equilibrium (Fi,0) — (R^e^^O) described by equations can be expressed using 
(xi,X2,ii,i2,i3) e according to L3.26\) - ^3^) . 
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Summarizing the above, the linearization of (|4.3p - l|4.4p about the inverted equilibrium (Fj, 0) is expressed 

as 

xi — mglixi — 0, (4-25) 
X2 - mgl2X2 — 0, (4.26) 
±3 = 0. (4.27) 

Note that the inverted equilibrium of the Lagrange-Poincare reduced equations has two negative eigenvalues, 
two positive eigenvalues and a zero eigenvalue. Thus, the inverted equilibrium (F^jO) is unstable and locally 
there exists a two dimensional stable manifold, a two dimensional unstable manifold and a one dimensional 
center manifold. 

Proposition 14. The inverted equilibrium (Fi,0) = (^J^jO); Lagrange-Poincare reduced 

dynamics of the 3D pendulum described by equations 114. 3\ l and |^.^[ ) is unstable. 

Remark 3. Note that combining Proposition [HI with Proposition [TO] immediately confirms the result 
that the inverted equilibrium manifold I of the full equations for the 3D pendulum given by (|3.ip - (|3.2p is 
unstable. 

We have analyzed the local stability properties of the hanging equilibrium and of the inverted equilibrium. 
We have not analyzed local stability properties of any other equilibrium solutions, but this analysis can easily 
be carried out using the methods that have been introduced. 

5. Lagrange— Routh Reduced 3D Pendulum Dynamics on TS^. In the previous sections we 
studied the full and the Lagrange-Poincare reduced equations of motion of the 3D pendulum. These involved 
the study of the dynamics of the 3D pendulum on TSO{3) and on TS0{3)/ , respectively, using {R,lu) 
and (F, lo) to express the equations of motion. In this section, we present Lagrange-Routh reduction of the 
3D pendulum, and we study the equations of motion that describe the evolution of (F,F) e TS^. 

5.1. Lagrange— Routh Reduction of the 3D pendulum. The key feature of Lagrange-Routh re- 
duction is reducing the configuration space into the quotient space induced by the symmetry action. The 
resulting equations of motion on the reduced space are described in terms of the Euler-Lagrange equation, 
but not with respect to the Lagrangian itself but with respect to the Routhian [15l [TU [12] . 

The 3D pendulum has a symmetry given by a rotation about the vertical axis. The symmetry action 
$e : 5*0(3) ^ 5*0(3) is given by 

<^g{R) = exp{9e3)R, 

for 6 Q and R G 5*0(3). It can be shown that the Lagrangian of the 3D pendulum is invariant under 
this symmetry action. Thus, the configuration space is reduced to the shape manifold 50(3)/5^ = 5*^, and 
the dynamics of the 3D pendulum is described in the tangent bundle TS^. This reduction procedure is 
interesting and challenging, since the projection 11 : 50(3) given by n(i?) = R^e^ together with the 

symmetry action has a nontrivial principal bundle structure. In other words, the angle of the rotation about 
the vertical axis is not a global cyclic variable. 

Here we present expressions for the Routhian and the reduced equations of motion. The detailed de- 
scription and development can be found in the Appendix. 

Proposition 15 ([15j). We identify the Lie algebra of with R. For {R,u!) £ TiiSO{3), the mo- 
mentum map J : TSO{5) K*, the locked inertia tensor I{R) ; M ^ M*, and the mechanical connection 
A : T5'0(3) — > R for the 3D pendulum are given as follows 

3{R,Lo) = elRJuj, (5.1) 
1{R) = ejRJR^es, (5.2) 
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The value of the momentum map fi = 3{R,lu) corresponds to the vertical component of the angular 
momentum. Noether's theorem states that the symmetry of the Lagrangian implies conservation of the 
corresponding momentum map. This is an alternative method of showing the invariant properties of the 3D 
pendulum, as opposed to the direct computation used in Section [3. II 

Based on the above expressions, Lagrange-Routh reduction is carried out to obtain the following result. 

Proposition 16. For a given value of the momentum map fi, the Routhian of the 3D pendulum is given 

by 

R'^ir^t) = i(r X r) • j(r X r) - ^{b^ + v^)ir ■ jt) + mgv ■ p, (5.4) 

where b = ^^^-w-jtt-^ , v — Y~Tr' '^'^'^ magnetic two form can be written as 



/3^(r X ,7,r X c) = - (fTjfyi [-(^ ■ J^MJ] + 2 ll^rfj r ■ (77 x c). (5.5) 

The Routhian satisfies the Euler- Lagrange equation, with the magnetic term, given by 

sf R''{T,t)dt= [ if,l3^{ST)dt. (5.6) 
Jq Jo 

This yields the reduced equation of motion on TS^ : 

f = -||r||2r + r X E, (5.7) 



wh 



ere 




(j(r X r) - bjT) X ((r x r) - bv) + v^jt x r - mgV x p - cr , (5.8) 



r • jr J ' T - JT ' r • jr 

Proof. Sec the Appendix. □ 

5.2. Lagrange Routh Reconstruction of the 3D pendulum. For a given value of the momentum 
map /i, let V{t) G 5^ be a curve in the reduced space S*^ satisfying the Euler-Lagrange equation for the 
reduced Routhian R^^ given by (|5.7p . The reconstruction procedure is to find the curve R{t) G 5*0(3) in the 
configuration manifold that satisfies n(i?(i)) = T{t) and i {R{t) , R{t)'^ R{t)) = fi. 

This is achieved in two steps. First, we choose any curve Rhor{t) G 50(3) in the configuration manifold 
such that its projection is equal to the reduced curve, i.e. n(i?hor(0) = ^(t). Now, the curve R{t) can be 
written as R{t) = $g(j)(i?hor(i)) for some 6{t) G 5^. We find a differential equation for 9{t) so that the value 
of the momentum map for the reconstructed curve is conserved. 

Proposition 17. Suppose that the integral curve of the Lagrange-Routh reduced equation \5. 7| j is 
given by {T{t),T{t)) G TS'^, and the value of the momentum map /i is known. The following procedure 
reconstructs the motion of the 3D pendulum to obtain {R(t),u){t)) G T5'0(3) such that \\{R{t)) — T{t) and 

J{R{t),L0{t))=fi. 

1. Horizontally lift T{t) to obtain _Rhor(i) by integrating the following equation with _Rhor(0) = i?(0). 

Rhorit) = Rhor{t)uJhor{t), (5.10) 

where 

whor(i) = Tit) X r{t) - bit)T{t). (5.11) 
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2. Determine 0dyn(i) £ 5*^ by the following equation. 

^<^-(^) = /r>r^^- (5.12) 

3. Reconstruct the curve in TS0{3). 

R{t) - $e,^„(t)(i?hor(0) - exp[0dy„(i)e3]-Rhor(i), (5.13) 

Uj{t)^UJi,or{t)+l^{t)T{t). (5.14) 

Proof. See the Appendix. □ 
This leads to the geometric phase formula that expresses the rotation angle about the vertical axis along 

a closed integral curve of the reduced equation. 

Proposition 18. Let T(t) be a closed curve in S^, i.e. T{0) ~ r(r) for some T. The geometric phase 

OgeoiT) G of the 3D pendulum is defined by the relationship R{T) = ^^(y)(i?(0)) when /x = 0. This can 

be written as 

« m /• 2||Jr(0f -tr[j] im-jm) 

where B is a surface in with Tit) as boundary. 

5.3. Invariants of the Lagrange Routh Reduced Model. In this section we find an invariant 
of the motion for the Lagrange-Routh reduced model of the 3D pendulum, namely the total energy of 
the system. Note that the Lagrange-Routh reduced equations of motion are derived by elimination of the 
conserved vertical component of the body-fixed angular momentum. In later sections, we make use of the 
constant energy surfaces to understand the dynamics of the 3D pendulum. 

Proposition 19. The total energy 

E ^]^{t xV + {v - b)TYj{t xT+{v- b)T) - mgp^V (5.16) 

is constant along motions of the Lagrange-Routh reduced equations for the 3D pendulum. 

Proof. Substituting the reconstruction equations for the angular velocity (|5.1ip . (|5.14p into the total 
energy expression (|4.6p . we obtain (|5.16p . The time derivative of the total energy is given by 

E={txV+{v- h)rfj{f xT +{v -b)T + {u - b)t) - mgp^t. 

Substituting the reduced equation of motion (j5.7p into the above equation and rearranging, we can show 
that E = 0. □ 

5.4. Equilibria of the Lagrange-Routh Reduced Model. The Lagrange-Routh reduced model 
can be considered as the Lagrange-Poincare reduced model where the angular velocity is projected onto 
TrS'^. Thus, the equilibria structure of the Lagrange-Routh reduced model is equivalent to the Lagrange- 
Poincare reduced model, but it is represented by conditions on the reduced attitude Fc and the value of the 
momentum map fi instead of (FejWe). 

More explicitly, we study the equilibria structure of the Lagrange-Routh reduced model using (|5.7p , and 
we show it is equivalent to the families of equilibria presented in Proposition [9] 

Proposition 20. Consider the Lagrange-Routh reduced model of the 3D pendulum given by ( [5. 7| j. The 
equilibria (FcO) £ TS^ of the Lagrange-Routh model are given for /i G M as follows. 

1. The hanging equilibrium: ^jj^iO^ /i = 0, 

2. The inverted equilibrium: M — 



3. Two relative equilibria: 



J-'pW ) ' VlU 



p 
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4- One dimensional relative equilibrium manifolds: 



-sgn(a)|| ^,0 /i = ±sgna) /- -^n^Juc, 



(5.18) 



where — {J — -^Isxs) ^ P G corresponding to a G Li, i G {1, 2, 3, 4} and 

A = (-oo,0)U(Jj,^), L2^{0,j-J, /:3^(J-,-L), c^^ij-^,j-^). 

Also, sgn(-) denotes the sign function. 
The families of equilibria given in 115. 18\) converge to the hanging equilibrium and the inverted equilibrium 
when a — » 0, and they converge to the third equilibria given in ^5.17^ , when a — > ±oo. If the vector from 
the pivot to the center of mass p, lies on a principal axis, i.e. p x a = for some i £ {1, 2, 3}, then 115. 18\) 
can be rewritten as {(ei,0), (—6^,0)} for any G R. Furthermore, there exist additional equilibria under the 
following assumptions on the moment of inertia matrix J and the vector from the pivot to the center of mass 
P = [Pi P2 PzY- 

5. Ji, J2 and J3 are distinct and pt = 0. Then there exist one- dimensional relative equilibrium mani- 
folds: 



foriG {1,2,3}„ where p, = (^,-^,^), p, ^ (^,^,^), p, ^ (^,^,^) and 

7 e R. 
6. Ji = J2 ^ J3. 

(a) If pi = P2 = 0. Then there exist two-dimensional relative equilibrium manifolds: 



13' 



where q — (7, 5, j^lj-^ ) '"^^ 7, (5 £ R. 
(b) If P3 — 0. Then there exist one- dimensional relative equilibrium manifolds: 



r 



^,0) ^i = ±J^r^Jr, 



where r = {j^^, jf^^l) and 7 £ M. 
Proof. Substituting To = into (|5.7p - (|5.9p . we obtain a condition for an equilibrium To for /i as 

To X [v'^JTc X Te - rngFe x p] = 0. 

This is equivalent to 

[v'^JT^ x Te - mgV^ x p] = fca JFo (5.20) 

for a constant ^2 £ M. Taking the dot product of this and Fc implies that = k2T^ JT^. Since T^JTc > 
as the moment of inertia matrix J is positive definite and Fc £ S'^, it follows that ^2 = 0. Thus, (|5.20p is 
equivalent to 

ly'^ JTc + mgp = kiTe (5.21) 

for a constant fci £ R. Note that this is equivalent to the equilibrium condition for the Lagrange-Poincare 
reduced model given by (I4.14p : for any solution (Fo,fc, /ci) of ()4.14p . we can choose p, such that = 
1/2 = (^r'-i ' ^jr )^ ' which gives a solution of (j5.2ip . and vice versa. Thus, the equilibria structure of the 
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Lagrange-Routh reduced model is equivalent to the equilibria of the Lagrange-Poincare reduced model. For 
an equilibrium (FcWo) of the Lagrange-Poincare reduced model, the value of the momentum map at the 
corresponding equilibrium of the Lagrange-Routh model is given by 

fi = HT^JT,) = c<jjre(r^ JPo) (5.22) 

Substituting this into the equilibria presented in Proposition [51 we obtain the equilibria of the Lagrange- 
Routh reduced model. □ 

5.5. Local Analysis of the Lagrange-Routh Reduced Model on TS^. We showed that in case 
l-j, = 0, the Routh reduced model of the 3D pendulum has two isolated equilibria, namely the hanging 
equilibrium and the inverted equilibrium. These equilibria correspond to the disjoint equilibrium manifolds 
of the full equations of the 3D pendulum. 

We next focused on these isolated equilibria of the Routh reduced equations. Using Proposition \20\ 
the stability properties of the equilibrium manifolds of the 3D pendulum can be deduced by studying the 
Routh reduced equilibria for the case /i = 0. Compared to the Lagrange-Poincare reduced model, the Routh 
reduction procedure results in a set of complicated equations that are a challenge to analyze. 

We now present local analyses of the Routh reduced model of the 3D pendulum near the hanging 
equilibrium and near the inverted equilibrium. 

Consider the equations (I3.22p - (|3.24p representing the linearization of the full equations of motion of 
the 3D pendulum at the hanging equilibrium. It was shown before that the Lagrange-Poincare reduced 
equations of motion can be written in terms of (xi, 2:2, cci, 12, is). As shown in Proposition llli this result 
follows from the fact that any perturbation in F € S"^ at can be expressed in terms of {xi, a;2) G K^. In 
a similar fashion, one obtains the following result. 

Proposition 21. The linearization of the Routh reduced attitude dynamics of the 3D pendulum, at 
the equilibrium (Fh,0) = (R^es^O), described by equation |5. 7| j can be expressed using (xi, 0:2, ii, ±2) G 
according to i3.2S\) and LS.23\) . 

Proof. In Proposition II H it was shown that the perturbations in F at Th can be described in terms of 
(xi, X2) e M^. The proof then simply follows by noting that the equations of motion of the Routh reduced 3D 
pendulum is described in terms of (F, F) S TS'^. Thus the linearization of the Routh reduced 3D pendulum 
model can be described using using (xi, X2, ii, ±2) S M'* according to (|3.22p and (|3.23p . □ 

The linearization of the Routh reduced attitude dynamics of the 3D pendulum, about the equilibrium 
(0,F^) is obtained from the linearized model of the full attitude dynamics by neglecting the dynamics 
corresponding to 0:3. Summarizing the above, the linearization of (15. 7|) about the hanging equilibrium 
(0, F/i) is expressed as 

Xi + mgliXi — 0, (5.23) 
X2+mgl2X2^0- (5.24) 

It is clear that due to the presence of imaginary eigenvalues, stability of the hanging equilibrium (F^,0) 
cannot be concluded. Therefore we next consider Lyapunov analysis. 

Proposition 22. The hanging equilibrium (F/i,0) = (^]j^jO) of the reduced dynamics of the 3D pendu- 
lum described by |5.7| ) is stable in the sense of Lyapunov. 

Proof. Consider the candidate Lyapunov function 

F(F, f ) = i(f X F + (i/ - 6)F) V(f X F + (t. - 5)F) + m5(||p|| - p^F). (5.25) 

Note that V^(Fh, 0) = and ^(F, F) > elsewhere. Furthermore, the derivative along a solution of (|4.3p and 
(|4.4p is given by 



y(F, f ) = (f X F + (i/ - b)Tyj{f xT+{v-b)T + {iy- b)t) - mgp'^t. 
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Substituting the reduced equation of motion (j5.7p into the above equation and rearranging, we can show 
that y(r,r) = O. Thus, the hanging cquihbrium of (|5.7[) is Lyapunov stable. □ 
Remark 4. Note that combining Proposition with Proposition [^D] immediately yields the result in 
Proposition O 

We next study the local properties of the Routh reduced equations of the 3D pendulum near the inverted 
equilibrium (ri,0). Consider the linearization of (|5.7p at an equilibrium (FijO) = (R^e^jO), where {Rc,0) is 
an equilibrium of the inverted equilibrium manifold I. A result similar to Proposition 1211 follows . 

Proposition 23. The linearization of the reduced attitude dynamics of the 3D pendulum, at the equilib- 
rium (FijO) = (i?Je3,0), described by |5. 7| j can be expressed using (zi, a;2, ii, G according to i3. 26]) 
and (3IT\) . 

Summarizing the above, the linearization of (j5.7p at the inverted equilibrium (ri,0) is expressed as 

Xi — mglixi — 0, (5.26) 
X2 — mgl2X2 — 0. (5-27) 

Note that the linearization of (j5.7p at the inverted equilibrium has two negative eigenvalues and two positive 
eigenvalues. Thus, the inverted equilibrium (ri,0) of the Routh reduced model is unstable and locally there 
exists a two dimensional stable manifold and a two dimensional unstable manifold. 

Proposition 24. The inverted equilibrium {Ti,0) = ("j^iOj the Routh reduced dynamics of the 

3D pendulum described by ( [5. 7p is unstable. 

Remark 5. Note that combining PropositionlMl with Proposition [501 immediately yields the result that 
the inverted equilibrium manifold I of the 3D pendulum given by (|3.ip ^ (|3.2p is unstable. 

5.6. Poincare Map on the Lagrange Routh Reduced Model. A Poincare map describes the 
evolution of intersection points of a trajectory with a transversal hypersurface of codimension one. Typically, 
one chooses a hyperplane, and considers a trajectory with initial conditions on the hyperplane. The points 
at which this trajectory returns to the hyperplane are then observed, which provides insight into the stability 
of periodic orbits or the global characteristics of the dynamics. 

The Lagrange-Routh reduced equations for the 3D pendulum on TS^ are a particularly suitable choice 
for analysis using a Poincare map, since it has dimension 4. Since the total energy given by (|5.16p is 
conserved, choosing a Poincare section on TS^ and restricting to an energy isosurface induces a Poincare 
map on a 2-dimensional submanifold of TS^. We define a Poincare section on TS^ for the Lagrange-Routh 
dynamics of the 3D pendulum given by (j5.7p as follows. 

7^ = |(r,r) G TS^ I ejt = 0, e|'(r X r) > 0, and E(r,t) = constant | . 

Suppose r G 7^ is given. The tangent space Tr^'^ is a plane that is tangential to and perpendicular to F. 
The first condition of the Poincare section, ejT — 0, determines a line in which the tangent vector F e TrS''^ 
should lie, and the constraint of the total energy conservation fixes the magnitude of the tangent vector in 
that line. Thus, the tangent vector is uniquely determined up to sign change. The second condition of the 
Poincare section resolves this ambiguity. It also excludes two reduced attitudes F = ±63 for which the first 
condition is trivial; ejF = for any F £ T^^S^ (jT-ggS*^. Thus, V can be equivalently identified as 

7^ = |f e 5^ I ejf = 0, 63 (F X f ) > 0, and E{T, f ) = constant | , 

where (F,f) satisfies (|5J| . 

This Poincare section in TS^ is well-defined in the sense that for each element, the corresponding tangent 
vector is uniquely determined. The attitude and the angular velocity in TSO{i) can be obtained by using 
the reconstruction procedure for the given value of the momentum map /i. 

Fig. 15.11 shows particular examples for this Poincare map. The pendulum body is chosen as an elliptic 
cylinder with properties of m = 1kg, J — diag[0. 13, 0.28, 0.17] kgm^, p — [0, 0,0.3] m. The initial condition 
are given by i?o = /sxs and ujo ~ c[l, 1, 1] rad/s, where the constant c is varied to give different total 
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(a) E = -2.65 (h) E = (c) E = 2.03 (d) E = 8.83 (e) E = 11.95 

Fig. 5.1. Poincare maps for 3D pendulum with varying total energy 

energy levels. The Lie group variational integrator introduced in |13j is used to compute the Poincare maps 
numerically. 

It is interesting to see the transition of the Poincare maps with varying total energy levels. The attitude 



dynamics of the 3D pendulum is periodic in Fig. 5.1(a) but it exhibits chaotic behavior with increased 



energy level in Fig. 5.1(b) and 5.1(c) If the total energy is increased further, the attitude dynamics becomes 



periodic again in Fig. 5.1(e) This demonstrates the highly-nonlinear, and perhaps chaotic, characteristics 



of the 3D pendulum dynamics. 

6. Conclusions. The 3D pendulum exhibits rich dynamics with nontrivial geometric structure; these 
dynamics are much richer and more complex than the dynamics of a ID planar pendulum or a 2D spherical 
pendulum. This paper has demonstrated that the methods of geometric mechanics and the methods of 
nonlinear dynamics can be meshed to obtain insight into the complex dynamics of the 3D pendulum. 

We have introduced three different models of the 3D pendulum, including the full model defined on 
TSO{3), the Lagrange-Poincare reduced model on TS0{3)/ obtained by identifying configurations in the 
same group orbit, and the Lagrange-Routh reduced model on TS^ where one additionally utilizes the fact 
that the dynamics evolves on a momentum level set. Relationships between the various representations are 
discussed in the context of conservation properties, equilibria and their stability properties, and invariant 
manifolds. 

In addition, we illustrate that the use of the Lagrange-Routh reduced equations of motion, together with 
the energy conservation properties, allow the construction of a Poincare map that can be readily visualized, 
thereby providing a graphical tool for obtaining insight into the rich nonlinear dynamical properties of the 
3D pendulum. 

Appendix. 

In this appendix, we summarize Lagrange-Routh reduction and reconstruction procedures for the 3D 
pendulum. 

A.l. Reduction. A description of Lagrange-Routh reduction can be found in [T5] including expressions 
for the mechanical connection and the Routhian of the 3D pendulum given by (|5.3p and (|5.4p . respectively. 
Here we derive the reduced equation of motion (|5.7p using the Euler-Lagrange equation for the given Routhian 

(lEl. 

Variation of Routhian. The Routhian satisfies the variational Euler-Lagrange equation with the magnetic 
term given by (|5.6p . We use constrained variations of P G 5^: 

ST^Txr], (A.l) 

dr = rxr] + rxfi. (A.2) 



Here we assume that 77 • P = 0, since the component of rj parallel to P has no effect on ST. These expressions 
are essential for developing the reduced equation of motion. 
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Using l|A.l[) . (|A.2[) . and the properties F • f = 0, F ■ 77 = 0, the variation of the Routhian is given by 

Ji?^ = ,) . J(f X F - 6F) - r; ■ F X -fx J(f x F) + {b^ + j/^) JF - 6J(f x F) + b{t x JF) + mgp . (A.3) 

Magnetic 2-form. From the given mechanical connection A and a value of the momentum map /i G K* , 
define a 1-form A^, on TSO{3) by 

A,{R) ■ {R,u;) ^ = t^^^^^ 

The magnetic 2-form /3p in (jS.Sp is the exterior derivative of Afi, which can be obtained by using the identity 
d^^(X,y) = X[A^{Y)]-Y[A^{X)]-A^j,{[X,Y]) forX = Rfi,Y = EC G TrS0{3). Suppose that f = Fxw. 
Since T ■ (uj x rj) = rj ■ (T x lu) = rj ■ T, the interior product of the magnetic 2-form is given by 



itf3^i6T) =p^{Txu;,rxr)) = i^{ tr[J] 



F • JF 



F • TJ, 



(A.4) 



where v — y~Jt- 



Euler- Lagrange equation with magnetic 2-form. Substituting (jA.3[) and (jA.4[) into (j5.6p . and integrating 
by parts, the Euler-Lagrange equation for the reduced Routhian ()5.4p is written as 



77 



J{t X T - bV - bT) + T X X + ct 



dt = 0, 



where 

X = -f X J(f X F) + (6^ + iy2^ jr - bJ{t X F) + 6(f X JF) + mgp, 
and c is given by (|5.9p . Since (|A.5p is satisfied for all ?/ with F • 77 = 0, we obtain 

J(f X F - fef - 6F) + F X a: + cf = AF, 



(A.5) 



(A.6) 



(A.7) 



for A G M. This is the reduced equation of motion. However, this equation has an ambiguity since the value 
of A is unknown; this equation is implicit for F since the term b is expressed in terms of F. The next step is 
to determine expressions for A and b using the definition of b and some vector identities. 

We first find an expression for A in terms of F, F. Taking the dot product of (jA.7P with F, we obtain 



F • J(f X F - 6f - &F) = A. 



(A.8) 



From the definition of 6, we can show the following identity: F • J(f x F — bV) = 0. Differentiating this with 
time, and substituting into (jA.Sp . we find an expression for A in terms of F, F as 



A = -F • J(F X F - bV). 



(A.9) 



Substituting (|A.9p into (jA.7p . and taking the dot product of the result with F, we obtain an expression for 
b in terms of F, F as 



6 = F • J-^ |f X X + cf + (f • J(f X F - 6F))f| 



(A.IO) 



Substituting (|XTOl) into (|X7|, and using the vector identity y- (F-r)F = (F-F)r- (F-r)F = -F x (F x F) 
for any K G M'^, we obtain the following form for the reduced equation of motion 



FxF-6F-Fx 



F X J-i |f X a: + cf + (f • J(f X F - 6F))f| 



= 0. 
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Reduced equation of motion. This equation has no ambiguity. Now, we simphfy this equation. The 
above expression is equivalent to the foUowing equation 



r X 



rxr-6r-rx 



r X |r X X + cr + (r • j(r X r - 5r))r| 



= 0. 



Since TV 



r|p. the first term is given by 

r X (f X r) = (r ■ r)f - (r • f )r 



iriPr. 



Using the property T x (T x (F x F)) = -(F • r)r x y -T x F for F e M^, the third term of the above 
equation can be simphfied. Substituting (|A.6[) and rearranging, the reduced equation of motion for the 3D 
pendulum is given by 



where S = 5r + J" 



r = -||r||^r + r x e, 

( J(f X r) - bJT) X ((f x r) - bT) + 1/2 jr x r - mgT X p-cT 



(A.ll) 



A. 2. Reconstruction. For a given integral curve of the reduced equation (r(i),r(t)) G TS^, we find 
a curve R{t) G SO{3) that is projected into the reduced curve, i.e. U{R{t) — r{t). The reconstructed curve 
can be written as R{t) = <I>g(()(i?(t)) for some 9{t) G S^. The conservation of the momentum map yields the 
following reconstruction equation [TS] . 

o{t)-^e{t) = r\R{t))^i - A{R{t)). 

The particular choice of R{t), the horizontal lift given by ()5.10p . simplifies the above equation, since the 
horizontal part of the tangent vector is annihilated by the mechanical connection. Further, since the group 
is abelian, the solution reduces to a quadrature such as (|5.12p . The reconstructed curve is given by (|5.13p . 
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